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Abstract 

We consider recurrent contagious processes on a time-varying network. As a control pro¬ 
cedure to mitigate the epidemic, we propose an adaptive rewiring mechanism for temporary 
isolation of infected nodes upon their detection. As a case study, we investigate the network of 
pig trade in Germany. Based on extensive numerical simulations for a wide range of parameters, 
we demonstrate that the adaptation mechanism leads to a significant extension of the parameter 
range, for which most of the index nodes (origins of the epidemic) lead to vanishing epidemics. 
Furthermore the performance of adaptation is very heterogeneous with respect to the index node. 
We quantify the success of the proposed adaptation scheme in dependence on the infectious pe¬ 
riod and detection times. To support our findings we propose a mean-field analytical description 
of the problem. 
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I. INTRODUCTION 


Recently the availability of data on host mobility and contact patterns of high resolu¬ 
tion offers many opportunities for the design of new tools and approaches for modeling 
and control of epidemic spread |[I]-l5]]. To exploit these versatile concepts of complex 
networks research, hosts or their spatial aggregations are considered as nodes and host 
contacts or their relocations as edges. Very frequently the edges are not static, but 
changing with time. If dynamical processes on networks possess a characteristic time 
scale much faster than the time scale of the changing edges, a static, quenched approx¬ 
imation of the topology may be sufficient. However if the time scale of the process is 
comparable with the time scale of the network change, a more sophisticated concept 
of time-varying or temporal networks is required, because static approximation might 
violate the causality principle ll6l-[T4ll. Although for static networks a variety of surveil¬ 
lance and control approaches were proposed based on various network measures such as 
node degree, betweenness centrality etc. IfTSlI . control concepts for temporal networks 
are still missing. The previous studies on this topic were devoted mostly to targeted vac¬ 
cination policies ifTfil - llTIl and general controllability questions Il22ll2^ . Furthermore, an 
adaptation of edges was proposed to mitigate the spread in static networks Il24l - l28l and 
by random rewiring, a co-evolution of the network and the spreading was implemented 
to avoid infected nodes. However, there have been no studies combining the adaptation 
approach of epidemic control with intrinsic temporal changes of the underlying network 
structure. 

In this paper, we propose an adaptive non-targeted control mechanism for spreading 
processes on temporal networks and assess its effectiveness. We consider a determin¬ 
istic recurrent contagious dynamics similar to a susceptible-infected-susceptible (SIS) 
model. In our model, susceptible nodes, after contact with an infected node, become in¬ 
fected and after a fixed time k they become again susceptible to the disease. We assume 
that the nodes are screened for infection, but the information whether or not a node is 
infected is available only after some detection time 6 . One can interpret the detection 
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time 5 as a time required to reliably detect the disease (also called window period |l29||) 
or a time the disease needs to manifest itself (incubation time) and be diagnosed. After 
the detection we apply adaptation rules, rewiring our system in a way to avoid edges 
emanating from the detected infected nodes. This results effectively in a temporary 
quarantine of the infected nodes.. 

We pursue the question, if the interplay between the intrinsic dynamics of a temporal 
network and adaptation rules leads to a substantial improvement of the disease mitiga¬ 
tion. In contrast to a range of studies considering targeted intervention measures, we 
apply control measures to all infected nodes, after they are detected as those. Our ap¬ 
proach can be easily supplemented by targeted interventions as well, where only some 
fraction of specific nodes is controled. 

This paper is structured as follows: First, we introduce the empirical dataset and 
outline our approach to adaptive epidemic control. Then, we present results on the 
mitigation strategy and discuss their implications. Finally, we summarize our findings 
and provide an outlook on further research directions. 


II. METHODS AND DATASET 

The empirical temporal network investigated in our study is extracted from the 
database on pig trade in Germany HI-Tier (See also Data Accessibility Section). We 
use an excerpt from this animal (pig) trade network with 15,569 agricultural premises 
(nodes) over a period of observation of 2 years (daily resolution), which consists of 
748,430 trade events (links). All premises have been anonymized for the study. On 
average day (except Sunday) the network contains 1220 nodes with 1141 edges as de¬ 
picted in Fig. 2(a) (Supplementary Material). We observe a non-uniform activity with 
respect to the day of the week (see Supplementary Material, Fig. 1): on average, we 
find 1346 nodes with 1258 edges on a working day and 932 or 43 nodes with 844 or 
27 edges on Saturday and Sunday, respectively. This reduced activity on the weekend 
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is clearly visible in Fig. 2(a) (Supplementary Material). The trade flow is directed: 
from source nodes (piglet producers) to sink nodes (slaughter houses). There are 291 
sinks, which only have in-coming links, but no out-going links, and 5,504 sources with 
only in-coming links. For fuhrer details on basic network characteristics, see Tab. 1 in 
Supplementary Material. 

The considered network also possesses a strong heterogeneity with respect to the size 
of out-components as shown in Fig 2(b) (Supplementary Material). The out-component 
of a node is defined as the number of nodes that could be infected in a worst-case SI 
(susceptible-infected) epidemic scenario with an outbreak originating from that particu¬ 
lar node upon its first occurrence. For this worst case, we assume the infinite infectious 
period and consider an SI epidemic following the directed, temporal links during the 
whole observation time. The distribution of the out-components peaks around 6,000 
nodes, which corresponds to 40% of the network. We also find many cases, where 
an outbreak immediately stops and the length of the respective epidemic path is short. 
Among these nodes are the above-mentioned sink nodes. 

Therefore, it can be expected that the prevalence, i.e. the total number of infected 
nodes at a given point in time, strongly depends on the outbreak origin, and a surveil¬ 
lance strategy that randomly selects nodes for screening will not be effective. In ad¬ 
dition, for finite infectious periods, the day of first infection is important as has been 
shown in Ref. IfTOl . 

In our simulations on this real-world temporal network, we consider a determinis¬ 
tic recurrent epidemic of an SIS type. Note that SIS epidemics on temporal networks 
without adaptation have been investigated in Ref. Ool . The spreading process is deter¬ 
ministic in the following sense: every time a trade event from an infected to a susceptible 
premise takes place, the susceptible one becomes infected with probability 1. In other 
words, we consider diseases with high infectiousness and thus a worst-case scenario 
for a spreading process. More precisely, nodes can be susceptible (S), infected unde¬ 
tected (/), or infected detected (D) as depicted in Fig[T| After detection time 6, infected 
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Figure 1. (Color online) Infection model and adaptation mechanism. Nodes are arranged ver¬ 
tically with links between them for each day. Nodes can be in three states: susceptible, S 
(empty circles), undetected infected, / (yellow), and detected infected, D (red). After 6 days, 
infected nodes are detected and for f G [C + S, + k] - with C being the time of infection - 
all out-going links from the detected nodes (red) are randomly rewired to other susceptible or 
undetected nodes as starting point, e.g. the link a is destroyed and instead the link b is created. 
Note that a newly created links might spread the disease as well. E.g. on day 5 the green arrow 
might point from the undetected infected node E to the node C (link c) spreading the disease. 
After K. days, infected nodes become susceptible again. Parameters: 5 = 2 d, = 4 d. 


nodes are detected and according to our control strategy, all of their out-going links are 
rewired to start at other susceptible or infected, but not yet detected nodes chosen at 
random. Thus, we isolate the detected infectious nodes or place them under quarantine. 
The total period of infection is denoted hy k. > 6. Note that nodes, which are infected, 
but not yet detected, take part in the rewiring and thus, eventually increase the risk of 
receiving nodes. See, for instance, the newly formed (green) link from node E to C at 
f = 5 d in Fig. The proposed adaptation scheme can be easily implemented for SIR 
models, where the infected nodes become immune, that is, recovered, after the infection 
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period. This would, however, considerably reduce the pool of available nodes over the 
course of the available observation period. 


III. RESULTS AND DISCUSSIONS 


In this section, we discuss the main findings based on the model and data described 
above. A typical evolution of the number of infected nodes (prevalence) is shown in 
Fig. 1^ for an arbitrary starting node and the infectious period k = 45 d. For a free- 
running sustained disease, i.e. without network adaptation or other measures of mitiga¬ 
tion (k = 5), one can observe that the prevalence fluctuates around a constant endemic 
level after a short transient period as shown by the blue curve. The effect of network 
adaptation manifests itself either in a reduction of this endemic prevalence level (orange 
curve in Fig. where network adaptation takes place after 5 = 25 d) or termination of 
an outbreak as shown by the red curve for 5 = 13 d. This shows that the considered 
adaptation scheme can substantially limit the spreading potential of an outbreak. This 
is also confirmed by the mean-field approximation (Section |IV]). 

In general, different nodes as origins of infection lead to different outcomes, some 
of them lead to epidemics and some do not. Note that we always start our simulation 
with just one origin (or index) node infected. In order to evaluate the influence of 
every possible node as an origin of infection, we scan the whole network by separately 
considering each node as the origin of an outbreak upon the node’s first appearance 
in the dataset. Technically we define epidemics as persistent if after 700 days (at the 
end of the maximal observation time) there still exists a non-zero prevalence in the 
system. Therefore we define the endemic fraction or the probability of an epidemic to 
be sustained as the fraction of index nodes leading to persistent epidemics 

.^endemic 


£ = 


N, 


( 1 ) 


total 


where Atotai is the total number of origin nodes and and Aendemic is the number of origin 
nodes leading to persistent epidemics. Without adaptation, £ increases with infectious 
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Figure 2. (Color online) Typical time course of an epidemic. Prevalence (number of infected 


nodes, daily resolution) for a fixed infectious period k = 30 d and different detection times 


6. Whereas unadapted epidemic with S = k. (blue curve), adapted epidemics with S = 37 d 


(purple), 5 = 25 d (orange) are persistent, adapted epidemics with 5 = 13 d (red) and 5 = 1 d 


(green) die out before the end of the observation. 




Figure 3. (Color online) Endemic fraction 6: in dependence on the infectious time k without 
control adaptive rewiring. 


period k. (Fig. [^. Only for very small infectious periods (/^ < 3 d) there is an almost 
vanishing endemic fraction due to the low frequency of network contacts. Values larger 
than 10% are found for infectious periods larger than 5 d. The highest endemic fraction 
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Figure 4. (Color online) Dependence of the endemic fraction £ on the infectious period k and 
detection time 5. In the region below the red line (5 ~ 0.14/v) the persistence fraction is less 
than 10%. 

in our system is around 70% for k = 200 d. The considered k- values fall in the biolog¬ 
ically plausible range including bacterial diseases such as hemorrhagic diarrhea caused 
by E. Coli with animals carrying the disease up to 2 months IIbTII . 

Besides the infectious period n the adaptation introduces an additional time scale 
— the detection time S. The dependence of the endemic fraction e in (k, 5)-parameter 
space is presented in Fig. The uncontrolled case of Fig. [^can be retrieved for S = k. 
We find that the (almost) disease-free region of small e becomes significantly larger than 
in the uncontrolled case, where the vanishing prevalence was observed for small k < 3. 
Due to the adaptation, the persistence fraction remains less than 10% for parameter 
values below the red line 5 0.14k -|- const. 

The behavior of the persistence fraction is also shown in Figs. 7(a) and (b) (Sup¬ 
plementary Material), which shows the persistence fraction e for fixed 6 or k, i. e., 
horizontal or vertical sections through Fig. respectively. In Fig. 7(a) (Supplemen¬ 
tary Material) we observe that in most cases for a fixed detection time the persistence 
fraction, starting from high values, increases intitially to eventually decrease with the 
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infectious period k. In Fig. 7(b) we see that for fixed infectious period k the endemic 
fraction monotonically increases with the detection time 5. 

For recurrent epidemics like an SIS process, the impact of the mitigation strategy by 
rewiring strongly depends on the node, from which the outbreaks originates, i.e. the 
index node. This heterogeneity of the prevalence is the direct implication of the het¬ 
erogeneity in the out-components (cf. Fig 2, Supplementary Material). The existence 
of clusters of nodes with similar prevalence values and adaptation properties is consis¬ 
tent with the existence of clusters of nodes with similar invasion routes as reported in 
Ref. ifT^ . There are subsets of index nodes, belonging to the same cluster, which are 
most responsive to adaptation. This knowledge can be exploited to target those nodes 
with high priority, if resources for disease control are sparse. To quantify the effect of 
adaptation on prevalence reduction for persistent epidemics, we define the efficacy as 


7 = 


(4 


r) 

ooJ 


( 2 ) 


where loo and denote the prevalence in the unadapted and adapted cases, respec¬ 
tively. Figure 5 (Supplementary Material) shows the influence of the index nodes and 
their impact on the success of the control measured in terms of the efficacy 7 for a fixed 
infectious period k = 28 d and different detection times S. For small S, we find the 
efficacy peaked at 90%. As the detection time becomes larger, the control becomes less 
effective indicated by smaller prevalence reduction. The mean efficacy averaged over all 
nodes as starting points is presented in Fig 8 in Supplementary Material. To character¬ 
ize the heterogeneity of the distribution of the 7 values shown in Fig. 5 (Supplementary 
Material), we also compute its entropy. Its dependence on the infectious period and 
detection time is presented in Fig. 9 (Supplementary Material). 
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IV. MEAN-FIELD DESCRIPTION 


We approximate the deterministic dynamics on a temporal network described in Sec¬ 
tion |^See Fig. 1) with the following system of stochastic reactions 


S +I ^21 
S + D ^21 



The first reaction describes the usual infection of a susceptible by an infective with the 
infection rate a which in the deterministic case could correspond to the average daily 
out-degree of a node a ~ {Kut)- The second reaction corresponds to the effective 
force of infection due to the occasional rewiring to infected nodes (See Fig. 1, link c) 
which reads a'SD with the effective infection rate a' = al/{I + S). The third reaction 
represent the detection with the rate v ~ 1/5. And the last one the recovery with the 
rate fi ~ l/(«^ — Using the fractions j = I/N, s = S/N, and z = D/N, the 
corresponding set of differential equations reads 



dz 




Note that the total number of nodes is conserved j + s + z = 1. Thus we could eliminate 
the third equation. Usually, to find the deterministic threshold for a disease outbreak, 
stability of the disease-free fix point (j = 0, s = 1) is considered. The eigenvalues are 


Al — —/i — -- 

5 — K 


and 


X 2 = a — 1/ = a — 1/5. 
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Because S < k, the first eigenvalue is always negative. From the outbreak condition 
A 2 > 0, we have aS > 1, i.e. is the standard threshold condition, without any depen¬ 
dence on the infectious period k. It seems that this mean-field approximation does not 
reproduce the observed threshold S ~ const x k -I- const correctly (see Fig. 4). The 
endemic prevalence is given by 



Note that in the limit /i —?• 00 or equivalently S = k, corresponding to the unadapted 
case we recover the well-known result for an SIS model: j = {a — v)/a. Thus the 
efficacy Q reads 



The endemic values of s* and z* variables are given by the expressions 



V. CONCLUSIONS 

We have investigated disease control based on an adaptive rewiring strategy of a 
time-varying network to mitigate the effect of a recurrent deterministic epidemic. This 
control measure relies on isolating infectious nodes and thus is different from most of 
control approaches proposed for temporal networks ifTbl - tT^I^ . We have considered a 
SIS-type dynamics on the nodes and introduced a detection time, after which links can 
be rewired to isolate infectious nodes. As an exemplary temporal contact network with 
real-world application, we analyzed an animal trade network, where each trading event 
corresponds to a contact between two agricultural premises. The network of farms can 
be seen as a contact network with nodes in a susceptible or an infected state. 

We have found that for recurrent epidemics, the starting point of an outbreak is very 
important for the course of the epidemics: it either dies out or becomes endemic with 


11 





different prevalence levels. This happens due to the heterogeneity of the subset of the 
network reachable from the specific first (index) node. Accordingly, we have found 
that the impact of a mitigation strategy by network adaptation is similarly variable. 
The region of disease parameters, where most of the index nodes lead to vanishing 
epidemics, can be substantially extended using the proposed adaptive rewiring strategy. 
To effectively control the epidemic, the detection times should be less than 10 days. 
Moreover, there is a range of detection time values between 7 and 10 days, which lead 
to especially effective mitigation of epidemics with an infectious period around 30 days. 
This might be due to the interplay of the internal time scales of the system. We have 
shown that the success of an adaptation depends also on the parameters of the epidemics 
and, for instance, saturates for very long infectious periods. 

In the presented work, we have provided a proof of concept and reported on the 
effect of modification of the contact network. The model can be further detailed and 
extended following a metapopulational approach, which takes into account the number 
of animals traded or present in the premises, heterogeneity of parameters, as well as 
stochastic effects along the lines of, for instance. Ref. ll3Tl - lT7]l . This, however, is beyond 
the scope of this study, but a promising topic for the future. 
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SUPPLEMENTARY MATERIAL 


The network of animal trade is visualized in Fig. The network exhibits a tree¬ 
like structure and can be interpreted as a hierarchical supply-chain network. This is 
of particular interest for the spread of diseases, because outbreaks will have different 
impact depending on where they first occur. 


Table I. Basic network properties of en exerpt of the time-aggregated German pig-trade network. 
Properties denoted by asterisk are for the largest connected component of the network considered 
as undirected. 


Property Value 

Number of nodes 15,569 

Number of edges 748,430 

Average daily number of nodes (except Sundays) 1220 
Average daily number of edges (except Sundays) 1141 
Average number of nodes (working days) 1346 

Average number of edges (working days) 1258 

Average number of nodes (Saturdays) 932 

Average number of edges (Saturdays) 844 

Average number of nodes (Sundays) 43 

Average number of edges (Sundays) 27 

Diameter* 13 

Average shortest path length* 4.081 

Average clustering coefficient* 0.1779 


In addition to Tab. basic network characteristics are presented in Figs. and 
Distribution of both the node degrees and the node activity aggregated over the whole 
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Figure 1. Visualization of the aggregated network of trading contacts, appearing at least 50 times 
in the dataset. Here the links are considered as undirected. 
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Figure 2. (a) Histogram of the number of active edges (blue) and nodes (green) (daily resolu¬ 
tion). (b) Distribution of the size of out-components for all nodes. The inset shows the same on 
a semi-log scale. 


17 











































































period of observation is very broad as shown in Fig.|^ 
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Figure 3. Histogram of node degree (blue: in-degree, green: out-degree) and activity (blue: 
in-coming activity, green: out-going activity) of the time-aggregated graph. 



Day of week 


LO 

c 


CD 

E 


1500 


1000 


500 


0 1——--u—--U— 

Mon Tue Wed Thu Fra Sat Sun 


Day of week 


Figure 4. (Average number of links and nodes on different week days. 

Figure [^depicts the average number of nodes and links resolved for each day of the 
week. Mondays show the highest numbers and other working days a similar level. On 
Saturdays, the network is less dense and the numbers of nodes and links are the lowest 
for Sundays. 

Figure shows the influence of the index nodes and their impact on the success of 
the control measured in terms of the prevalence reduction 7 ([T])for a fixed infectious 
period k = 32 d and different detection times S. For small S, we find a prevalence 
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Figure 5. Histogram of the prevalence reduction 7 for infectious period = 32 d and different 
detection times 6 as indicated in the panels. Each node in the network is chosen as a starting 
node in a separate simulation with initial infection upon its first appearance in the dataset. 
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Figure 6. Dependence of the endemic fraction e: left on the infection period k for different 
detection times S and right on the infectious period S for different detection times k. 


reduction peaked at 90%. As the detection time becomes larger, the control becomes 
less effective indicated by smaller prevalence reduction. The mean prevalence reduction 
averaged over all nodes as starting points is presented in Fig[^ 

Figs.j^a) and (b) shows the dependence of the endemic fraction £ ([T]) on k (for fixed 
6) and on 5 (for fixed k) correspondingly. 

Figure depicts the efficacy or prevalence reduction for sustained epidemics 7 Q 
in dependence on k and S. This quantity shows how useful the control is, even if we 
cannot totally stop the disease, which became endemic, and to what extent we could 
lower the prevalence. There is a fast decrease in the efficacy of the control measures as 
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Figure 7. Prevalence reduction £: in dependence on both infectious period and detection time. 
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Figure 8. Efficacy or relative prevalence reduction 7 averaged over all possible index nodes for 
sustained epidemics in dependence on k (with fixed 5, left panel) and on S (with fixed k, right 
panel). 

the detection takes longer. 

The efficacy or prevalence reduction 7 is highly dependent on the index node, and 
to characterize the heterogeneity we depict in Figure the entropy of the efficacy or 
prevalence reduction 7 defined as 

Sj = - J^7r(7i)log[7r(7i)], 

i 

where the index i enumerates all index nodes, in dependence on k and 5. 71 ( 7 *) is 
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Figure 9. Entropy in dependence on k and <5, which characterizes the heterogeneity in the distri¬ 
bution of efficacy 7 for different index nodes. 

the probability for an index node i to have the efficacy 7 *. Entropy in dependence 
on 5 (Fig. 1^ left panel) possesses clearly two flat levels, decreasing from a high to a 
low one with increasing S. Thus with early detection, different index nodes lead to 
very heterogeneous prevalence levels. Later detection leads to epidemics with similar 
efficacy levels. The dependence of the entropy on the infectious period k (Fig. right 
panel) exhibits maxima for the intermediate values of k. The largest heterogeneity in 
the efficacy for intermediate infectious periods might be due to the interplay of internal 
scales of the temporal network. 
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